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We study correlations between center vortices and the low-lying eigenmodes of the Dirac operator, in both 
the overlap and asqtad formulations. In particular we address a puzzle raised some years ago by Gattnar et al. 
[Nucl. Phys. B 716, 105 (2005)], who noted that the low-lying Dirac eigenmodes required for chiral symmetry 
breaking do not appear to be present in center-projected configurations. We show that the low-lying modes are in 
fact present in the staggered (asqtad) formulation, but not in the overlap and "chirally improved" formulations, 
and suggest a reason for this difference. We also confirm and extend the results of Kovalenko et al. [Phys. 
Lett. B 648, 383 (2007)], showing that there is a correlation between center vortex locations, and the scalar 
density of low-lying Dirac eigenmodes derived from unprojected configurations. This correlation is strongest 
at points which are associated, in the vortex picture, with non-vanishing topological charge density, such as 
vortex intersection and "writhing" points. We present supporting evidence that the lowest Dirac eigenmodes, 
in both asqtad and overlap formulations, have their largest concentrations in point-like regions, rather than on 
submanifolds of higher dimensionality. 



I. INTRODUCTION 

There are two reasons — "theoretical" and "experimental" 
— to believe that center vortices are the dominant feature of 
the vacuum state of pure SU(N) gauge theories at large scales. 
The theoretical reason is the simple fact that the asymptotic 
string tension of static color charge sources depends only on 
the Af-ality of the color charges. This fact can be understood 
as arising from color charge screening by gluons, which is an 
explanation in terms of particle excitations. But there should 
also be a corresponding "field" explanation of A^-ality depen- 
dence, in terms of gauge field configurations which dominate 
the path integral at large scales, and the center vortex scenario 
is the only explanation of this type which is known. The "ex- 
perimental" reason for believing in vortex dominance is the 
wealth of numerical evidence in its favor, which was summa- 
rized a few years ago in Ref. [1]. A key feature of this ev- 
idence is the strong correlation between vortex location, de- 
termined by center gauge-fixing and projection methods, and 
gauge-invariant observables, such as action density and the 
phase of Wilson loops. 

Center vortices were originally introduced to explain con- 
finement, but a force strong enough to confine quarks is also 
generally expected to break chiral symmetry [2]. 1 According 
to the Banks-Casher analysis [3], chiral symmetry breaking 
{%SB) is necessarily associated with a finite density of near- 
zero eigenmodes of the chiral-invariant Dirac operator, so we 
would of course expect this to be true for the Dirac spectrum 
evaluated in an ensemble of center-projected lattice configu- 
rations, which are known to be confining. 

Several years ago, however, Gattnar et al. [4] reported a 



The converse, of course, is not true. It is possible to have chiral symmetry 
breaking, as in the Nambu-Jona-Lasinio model, without having confine- 
ment. 



puzzling result. These authors computed the low-lying eigen- 
value spectrum of a chirally-improved version of the Dirac 
operator due to Gattringer [5], which approximates Ginsparg- 
Wilson fermions. A dense set of near-zero eigenvalues was 
found for unmodified lattice configurations, as expected, and 
a large gap in the eigenvalue spectrum, centered at eigen- 
value X = 0, opened up in the spectrum when center vor- 
tices were removed. This gap was also expected, since the 
vortex-removed configurations are not confining, and it has 
been known for a long time that (xjn|/) — > in vortex-removed 
configurations [6]. The puzzle was that an even larger gap in 
the spectrum was found for center-projected configurations, 
which contain only thin vortex excitations, and which are 
confining. Now if %SB was present in center-projected (or 
"vortex-only") configurations, and not in vortex-removed con- 
figurations, we would conclude that vortices explain %SB as 
well as confinement. If, on the other hand, %SB occured in 
vortex-removed configurations, but was absent in vortex-only 
configurations, we would then have to conclude that vortices 
are not especially relevant to %SB (although we would then 
like to understand how confinement can coexist with unbroken 
chiral symmetry). But the finding that %SB, or, to be more pre- 
cise, near-zero modes, are absent both in vortex-removed and 
in vortex-only lattices was unexpected, and it poses a chal- 
lenge to interpretation. 

The question which comes to mind is whether the large 
gap found by Gattnar et al. is related to the way in which 
chiral symmetry is realized on the lattice. The Casher argu- 
ment [2] that confinement implies %SB is based on the usual 
SU (N/)l x SU (Nf)u symmetry of the continuum theory with 
massless fermions. However, the chirally-improved Dirac op- 
erator only approximates this symmetry for gauge-field con- 
figurations which vary smoothly at the lattice scale. Center- 
projected configurations are not even close to smooth; pla- 
quette variables make a sudden transition from the trivial cen- 
ter element outside the thin vortex, to a non-trivial center el- 
ement inside. The chirally-improved Dirac operator is not 
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necessarily chirally symmetric, even approximately, in such 
backgrounds. In the absence of a symmetry, there is no rea- 
son to expect spontaneous symmetry breaking. If this fact 
explains why there is a gap in the eigenvalue spectrum of the 
chirally-improved operator, then it is reasonable to also expect 
a gap in the spectrum of the overlap operator [7], when evalu- 
ated on center-projected configurations. Of course the overlap 
operator, in contrast to the chirally-improved operator, does 
have an exact global symmetry, but the symmetry transfor- 
mations are gauge-field dependent [8], and only approximate 
the SU (Nf)^ x SU (Nf)g chiral symmetry transformations of 
the continuum theory for configurations which vary slowly at 
the scale of the lattice spacing. While this smoothness con- 
dition is expected in the continuum limit, it is never the case 
for center-projected configurations, and the Casher argument 
relating confinement to %SB need not apply. 

On the other hand, the Lagrangian for staggered fermions 
(and their asqtad cousins [9]) is known to be invariant un- 
der a subgroup of the usual chiral symmetry, irrespective of 
the smoothness of the gauge-field background. If the puz- 
zling gap in the Dirac eigenvalue spectrum found by Gattnar 
et al. is a consequence of the roughness of center-projected 
lattices, then we might expect this gap to disappear in the 
spectrum of the staggered or asqtad Dirac operators. Indeed, 
there is already a relevant result in ref. [10], which reported 
that (\jn|/) > for staggered fermions on a center-projected 
lattice. Likewise, suppose we somehow "soften" the center- 
projection procedure to make the center-projected configu- 
rations smoother, and evaluate the spectrum of the overlap 
operator on these smoothed vortex configurations. Then, if 
the roughness of thin vortices is the problem for the overlap 
formulation, the eigenvalue gap should go away for suitably 
smoothed, but still confining, vortex configurations. 

In section 2 we will report our results for the spectrum 
of the overlap and asqtad Dirac operators, when evaluated 
on normal, vortex-only (i.e. center-projected), and vortex- 
removed lattices. 2 Those results support the view that center 
vortices alone can induce both confinement and chiral sym- 
metry breaking. We go on in section 3, following the earlier 
work by Kovalenko et al. [13], to report on other correlations 
between center-vortex location, and the density distribution 
of low-lying Dirac eigenmodes in overlap and asqtad formu- 
lations. These correlations are an important test of the pic- 
ture advocated by Engelhardt and Reinhardt [14], in which 
topological charge is concentrated at points where vortices ei- 
ther intersect, or twist about themselves ("writhe") in a cer- 
tain way, and Dirac zero modes are concentrated where the 
topological charge density is large. If topological charge is 
concentrated in point-like regions, as is the case in the vortex 
picture, and if zero (and near-zero) modes are concentrated 
in regions of high topological charge density, then one would 
expect that the eigenmode densities of low-lying eigenmodes 



would be peaked in point-like regions. In section 4 we pro- 
vide some supporting evidence for this type of concentration. 
Our conclusions are found in Section 5. In an Appendix we 
review details of the tadpole-improved Ltischer-Weisz gauge 
action, and report on some necessary checks of vortex location 
via center projection, in numerical simulations of this lattice 
action. 



II. THIN VORTICES AND NEAR-ZERO MODES 

We begin with some preliminary information. Through- 
out this article we work with lattices generated by lattice 
Monte Carlo simulation of the tadpole improved Liischer- 
Weisz pure-gauge action, mainly at coupling $lw = 3.3 (lat- 
tice spacing a = 0.15 fm) for the SU(2) gauge group. The 
locations of center vortices are identified as usual by mapping 
the SU (2) lattice to a Z2 lattice which contains, by definition, 
only thin vortex excitations. The mapping is carried out by 
fixing the lattice to the direct maximal center gauge, which 
is equivalent to Landau gauge in the adjoint representation, 
and which maximizes the squared trace of link variables. The 
gauge-fixing procedure is the over-relaxation method. 3 The 
mapping to link variables on the center-projected (or "vortex- 
only") lattice, for the SU (2) gauge group, is given by 



U^x) -> Z^x) = signTr U^x) 



(2.1) 



and the link variables U' on the vortex-removed lattice are 
defined as 



U^x)=Z fl (x)U,(x) 



(2.2) 



The claim is that the thin vortices of the center-projected 
lattice lie somewhere in the middle of thick center vortices 
(thickness w 1 fm) on the unprojected lattice, and that thick 
center vortices are responsible for the area-law falloff of large 
fundamental Wilson loops, as well as for the A^-ality depen- 
dence of higher-representation string tensions in SU(N). The 
justifications for these claims, obtained from Monte Carlo 
simulation of the Wilson action, were reviewed some years 
ago in Ref. [1]. Similar tests in the case of the Ltischer-Weisz 
action are reported in the Appendix. Given that thin vortices 
locate thick vortices, then in the "vortex-removed" lattice we 
are really inserting a thin vortex somewhere in the middle of 
a thick vortex; the effect is to cancel out the field of the thick 
vortex at large distances. 

In Fig. 1 we display the twenty lowest-lying complex con- 
jugate eigenvalue pairs of the overlap Dirac operator, on a 
16 4 lattice at Plw = 3.3. Results are displayed for eigen- 



2 Similar results for the overlap spectrum on vortex-removed lattices have 
been obtained previously by Gubarev et al. [11] and by Bornyakov et al. 
[12]. 



3 In this study we perform 150 gauge fixing iterations for five gauge copies, 
and for the best copy continue with the gauge fixing for a further 400 it- 
erations. It is likely, however, that the selection of the best copy out of 
five copies is inessential, as the results are hardly distinguishable from a 
random choice of copy, as noted in the Appendix. 
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FIG. 1: The first twenty overlap Dirac eigenvalue pairs on the Ginsparg-Wilson circle for a 16 lattice, periodic boundary conditions at 
$lw = 3.3. The center-projected configurations show a four-fold degeneracy. 
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FIG. 2: The first twenty overlap Dirac eigenvalue pairs on the Ginsparg-Wilson circle for a 16 lattice, antiperiodic boundary conditions at 
Plw = 3.3. The zero modes in vortex-removed configurations disappear. 



modes obtained on the original, center-projected and vortex- 
removed lattices. We observe the same phenomenon already 
reported by Gattnar et al. for the spectrum of the chirally- 
improved Dirac operator: a large gap has opened in both the 
vortex-removed and the vortex-only spectra. We note here 
that the gap in vortex-removed overlap spectrum was found 
previously, and discussed in some detail, by Bornyakov et al. 
[12], who simulated a tadpole-improved Symanzik lattice ac- 
tion. Looking more closely at the spectra, we see that there 
only appear to be five eigenvalue pairs (out of twenty) in the 
center-projected case, indicating a four-fold degeneracy when 
the overlap operator is applied to Z2 lattice configurations. 
This factor of four has the following origin: In the first place, 
when link variables are simply plus or minus the 2x2 iden- 



tity matrix, the two colors decouple, and we have a factor of 
two degeneracy. Secondly, whenever the link variables are 
real and the Dirac operator has the Wilson or overlap (but not 
staggered) form, the eigenvalue equation D\\t„ = X„\\t„ is in- 
variant under charge conjugation. Thus, if \|/„ is an eigenstate 
with eigenvalue X n , then C \|/* is also an eigenstate, with the 
same eigenvalue [15]. This gives another factor of two, result- 
ing in an overall four-fold degeneracy. In the vortex-removed 
case, one does observe four near-zero modes of each chiral- 
ity, but these can be interpreted as a remnant of the exact zero 
modes of the free theory that are associated with a periodic lat- 
tice. These near-zero modes of the vortex-removed lattice are 
irrelevant to %SB, and disappear when we impose antiperiodic 
boundary conditions, as shown in Fig. 2. 
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FIG. 3: The first twenty asqtad Dirac eigenvalue pairs from a 16 4 lattice, periodic boundary conditions at fiuy = 3.3. 
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FIG. 4: The first twenty asqtad Dirac eigenvalue pairs from a 16 4 lattice, antiperiodic boundary conditions at fiiw = 3.3. 



In the Introduction we speculated that the reason for the 
large gap in the vortex-only case was connected with the 
lack of smoothness of center-projected lattices. In that case 
the exact symmetry of the overlap operator is strongly field- 
dependent, and does not really approximate the chiral symme- 
try of the continuum theory. Staggered and asqtad fermions, 
on the other hand, preserve a subgroup of the usual continuum 
SU(N/)l x SU(N/)r symmetry, irrespective of the smooth- 
ness of the configuration, and by the Casher argument [2] one 
would expect this remaining symmetry to be spontaneously 
broken, on the lattice, by any ensemble of gauge configura- 
tions with the confinement property. Then, according to the 
Banks-Casher relation, there should not be any gap in the 
vortex-only eigenvalue spectrum. 

In Fig. 3 we show our results for the twenty lowest-lying 
eigenvalue pairs of the asqtad Dirac operator, evaluated for pe- 
riodic and antiperiodic boundary conditions respectively. We 
see that the eigenvalue gap in the overlap spectrum, found on 
vortex-only lattices, is erased when eigenvalues are computed 



for the asqtad Dirac operator. In fact, the near-zero eigen- 
value density in the vortex-only case appears to be substan- 
tially enhanced, as compared to the density for unmodified 
lattices. In the vortex-removed case there also appear to be 
some near-zero modes centered at A. = 0, which are separated 
by a gap from the higher modes. However, a count of the 
number of these modes reveals that there are eight doubly de- 
generate eigenvalues in the central band with Im(A-) > 0, and 
an equal number of complex conjugates with Im(A-) < 0, for 
a total of 32 eigenmodes. Now, the free-field Dirac operator 
for massless staggered fermions has exactly four zero modes 
for each of four "tastes", and this number must be multiplied 
by the number of colors (i.e. two for SU (2)), for a total of 32 
free-field zero modes. So it is reasonable to guess that the 32 
eigenmodes in the central band of the vortex-removed spec- 
trum are simply the would-be zero modes of the free staggered 
theory. To check this, we carry out the same eigenmode calcu- 
lation using antiperiodic boundary conditions in one direction, 
which is sufficient to remove the zero modes of the free theory. 
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FIG. 5: The first twenty overlap Dirac eigenvalue pairs from a single configuration on a 16 4 lattice, antiperiodic boundary conditions at 
Plw =3.3, for interpolated fields. 



The result is shown in Fig. 4, where it is seen that the boundary 
conditions have no noticeable effect on the eigenvalue distri- 
bution for the original and center-projected lattices, while the 
32 eigenvalues of the central band completely disappear in the 
vortex-removed case. This result confirms the conjecture that 
the central eigenvalues for vortex-removed configurations are 
simply a remnant of the free-field zero modes, and play no 
role in %SB. Thus, for the asqtad operator, we have found ex- 
actly what was expected prior to the results of Gattnar et al.: 
the vortex excitations of the vortex-only lattice carry not only 
the information about confinement, but are also responsible 
for %SB via the Banks-Casher relation. This result was antici- 
pated in ref. [10], which found a non-zero (\j/V|/) condensate on 
center-projected lattices. %SB disappears for vortex-removed 



lattices, as discovered long ago by de Forcrand and D'Elia, in 
a direct calculation of (vjA|/) [6]. 

If the overlap operator yields misleading results on center- 
projected lattices, because of the lack of smoothness of center- 
projected configurations, then perhaps the overlap operator 
would produce a more reasonable answer when applied to a 
smoother version of the center-projected lattice. We therefore 
consider the following procedure: Given that SU (2) group el- 
ements can be represented by unit 4-vectors a A „ where U = 
aoh + icikO k , let Q/j(x) denote the angle between the vector 
representing group element U^(x) in maximal center gauge, 
and the vector representing the SU (2) center element Z (U (x)/2, 
where Zjx) was defined in eq. (2.1). Center projection simply 
takes this angle to zero, at every link, but we may also con- 
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FIG. 6: Finite temperature and center projection. The first twenty asqtad Dirac eigenvalue pairs from 12 3 x Nj center-projected lattices at 
Plw = 3.5 and Nj = 2,4,6, 12 lattice spacings, with the usual antiperiodic boundary conditions in the time direction. 



sider partial projections in which 8 ;U (x) is everywhere reduced 
by some fixed percentage. These partial projections interpo- 
late between the unprojected lattice, in maximal center gauge, 
and the fully center-projected lattice. In Fig. 5 we show the 
low-lying eigenvalues for partial projections, with 8 ;U (x) re- 
duced by 25%, 50%, 70%, 75%, . . ., together with the unpro- 
jected (0%) and fully (100%) center-projected lattices. We see 
that there is no really obvious gap in the partially-projected 
lattices, even at 80% projection. This agrees with our conjec- 
ture that applying the overlap operator to a smoother version 
of the vortex-only vacuum would give a result consistent with 
%SB and the Banks-Casher relation. 

We conclude this section with a high-temperature result. 
Chiral symmetry is restored at high temperatures, and this 
fact should also hold for center-projected lattices. There- 
fore, at sufficiently high-T, a gap should open in the eigen- 
value spectrum. This can be seen in Fig. 6, where we display 
the low-lying eigenvalues at (3/^ = 3.5 for time-extensions 
Nt = 2,4,6,12 computed on center-projected lattices. We 
note that the theory is certainly in the deconfined phase at 
Nt = 6, where there is, however, no obvious gap in the eigen- 
value spectrum, so it may be that on the projected lattice the 
chiral transition occurs at a higher temperature than the de- 
confinement transition. This is also consistent with ref. [10], 
which found a nonvanishing \jn|/ condensate at a temperature 
somewhat above the Wilson action deconfinement tempera- 
ture. We should stress, however, that there is no reason that 
the chiral and deconfinement temperatures need coincide on 
the center-projected lattice. Confinement is a sufficient but not 
a necessary condition for %SB, and, while the center-projected 
lattice is expected to get the static quark potential about right 
asymptotically, this fact certainly does not hold true at inter- 
mediate scales, where the finite thickness of real vortices is 
crucial. If the static potential on the center-projected lattice 
is strong enough, %SB will be realized, and this symmetry 
breaking may persist somewhat above the actual deconfine- 



ment temperature. The point is that if we view Yang-Mills 
configurations as being in some sense factorizable into vor- 
tices x perturbative fluctuations at short distances, then the 
thickness and internal structure of vortices is important for 
certain non-perturbative phenomena (such as Casimir scaling, 
and perhaps the precise chiral transition point) which are sen- 
sitive to the static quark potential at intermediate scales. 4 



III. VORTEX SURFACES AND DIRAC EIGENMODE 
DENSITIES 

The breaking of axial U(l) symmetry is associated with 
topological charge density in the vacuum state, and it is 
clearly of interest to understand the sources of this topolog- 
ical charge. There is a vast literature on this subject, and can- 
didate sources include instantons, calorons, and intertwined 
3-manifolds, with lattice-scale separation, of opposite topo- 
logical charge density [17]. It has also been suggested, by En- 
gelhardt and Reinhardt [14], that topological charge density 
is concentrated in certain regions of center vortices, where the 
vortices intersect or "writhe" (twist about) in some way. Since 
it is generally expected that zero modes of the Dirac opera- 
tor tend to concentrate in regions of large topological charge 
density, a correlation between the locations of vortex inter- 
section/writhing points, and the density p\{x) = \\yx(x)\ 2 of 
eigenmodes of the Dirac operator D, where D\\ix = X\\ix with 
X = (overlap formulation) or A. ~ (asqtad), would tend to 
support the Engelhardt-Reinhardt picture. A lack of corre- 
lation (or perhaps even an anti-correlation) would of course 
disfavor that picture. 



Cf. ref. [16] for a discussion of vortex thickness and Casimir scaling. 
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In our investigation we follow the approach of Kovalenko 
et al. in Ref. [13], who worked with eigenmodes of the over- 
lap operator, derived from lattices generated with the Wilson 
action. Kovalenko et al. proposed, as a measure of vortex- 
eigenmode correlation, the observable 

cM = Z P X x e H (V Px{x )-l) (3i) 
L Pi LxeH 1 

This choice of correlator requires some explanation. Center 
vortices on the full lattice are located, as explained previously, 
by center projection in maximal center gauge. Plaquettes on 
the projected lattice are either +1 or —1; plaquettes with the 
latter value are known as "P-plaquettes". However, the thin 
vortices of the center-projected configurations actually live on 
the dual lattice. In D = 4 dimensions each P-plaquette corre- 
sponds to a certain plaquette on the dual lattice, and these sets 
of plaquettes on the dual lattice form closed surfaces, namely, 
the thin center vortices. Each point on the vortex surface may 
belong to a certain number, N v > 3, of plaquettes on the vortex 
surface. If the surface is flat at the given point, then N v = 4. 
If the point is a corner of the surface, then N v = 3. When 
two flat vortex surfaces intersect at a point (or the same sur- 
face intersects itself), then N v = 8. It is also possible for the 
surface to twist (or "writhe") around a given point, in such a 
way that N v = 6. These writhing points are best visualized 
by consulting the illustrations in Ref. [14]. N v can take on 
other values as well, but N v = 1 or 2 is impossible for a closed 
vortex manifold. N v = holds for all points which do not 
belong to a vortex surface. For thin vortices, intersection and 
writhing points are clearly points where the topological charge 
density is non-zero, and the Engelhardt-Reinhardt proposal is 
that topological charge is also concentrated on the unprojected 
lattice in the neighborhood of these locations. In the definition 
(3.1) of the vortex-eigenmode correlator, the first sum is over 
all points p, on the dual lattice which belong to N v plaquettes 
on the thin vortex surface. The second sum over x G H is a 
sum over the 16 points in a hypercube H on the original lat- 
tice, surrounding the given point on the dual lattice. V is 
the lattice volume. 

Cx(N v ) is a measure of the average fractional excess of the 
eigenmode density, p^, over its average value (p^) = 1 /V, 
in the neighborhood of a point on the thin vortex connected 
to N v vortex plaquettes. If this value is non-zero, then the 
question is whether that value is large enough to be signif- 
icant. That is not a simple question to answer, in view of 
the fact that the location of thin vortices varies considerably, 
from one Gribov copy to another, so we certainly cannot deter- 
mine the vortex intersections of thick vortices with any great 
accuracy. This imprecision necessarily lowers the measured 
valued of the correlator Cx(N v ). The best we can do, at this 
point, is to compare the values of C\(N V ), with Dirac eigen- 
modes computed on the full lattice, to corresponding values of 
Cx(N v ), with Dirac eigenmodes computed on center-projected 
lattices. In the latter case, the only possible source of topolog- 
ical charge is the thin vortices, the topological charge density 
is highly localized, and the vortex location is certain. 

In Fig. 7 we display the data for Cx(N v ) vs. N v computed 



for eigenmodes of the asqtad Dirac operator in the full and 
(for comparison) center-projected configurations. The lattices 
are the same as in the previous section, generated by Monte 
Carlo simulation of the Liischer-Weisz action at fiiw = 3.3. 
Each plot displays the results for the first eigenmode (lowest 
X) and the twentieth Dirac eigenmode; we see that for the asq- 
tad eigenmodes there is not much difference. We find that the 
values of C\(N V ) obtained from eigenmodes in the full config- 
urations are only about a factor of four smaller than the corre- 
sponding values in the center-projected configurations, which 
in our opinion is not such a great suppression considering (i) 
the uncertainties in vortex location in unprojected configura- 
tions, i.e. intersection points of thin vortices may not precisely 
correspond to the actual intersections of thick center vortices; 
and (ii) the fact that center vortices are thin, i.e. extremely lo- 
calized, in the center-projected configurations, and therefore 
one expects a far greater degree of localization in the cor- 
responding Dirac eigenmodes, computed on center-projected 
lattices. But surprisingly, apart from the overall factor of four 
or so, the Figs. 7a and 7b look much the same. The most 
important feature, in our opinion, is the fact that the corre- 
lator increases steadily with increasing N v , and therefore the 
Dirac eigenmode density is significantly enhanced in regions 
of large N v . This fact seems at least compatible with the 
general picture advanced by Engelhardt and Reinhardt. We 
note, however, that the enhancement seems to be equal for all 
the low-lying eigenmodes, and not just the "would-be" zero- 
modes of the asqtad operator. 

Figure 8 shows our corresponding results for Cx(N v ) vs. 
N v , this time computed for eigenmodes of the overlap oper- 
ator. Our results in Fig. 8a are consistent with the previous 
results reported by Kovalenko et al. for overlap eigenmodes in 
Ref. [13]. In the case of the overlap Dirac operator there are 
true zero modes, and the correlation of the densities of these 
eigenmodes with vortices is significantly larger than the cor- 
relations of the 20th mode. In the case of the overlap operator 
it is pointless to compute eigenmodes on the center-projected 
lattice, for reasons we have discussed in previous sections. If 
we do this anyway, then we actually find an anf/correlation be- 
tween the low-lying modes and the vortex surfaces, as shown 
in Fig. 8b. 



IV. DIMENSIONALITY OF DIRAC EIGENMODE 
CONCENTRATIONS 

The data of the previous section provides some degree of 
evidence that low-lying Dirac eigenmodes concentrate pref- 
erentially at regions on the center vortex surface where there 
are, e.g., vortex self-intersections, and/or some sort of vortex 
twisting, such that "connectivity" N v of the vortex at a site is 
larger than N v = 4. Since surfaces intersect at points in four 
dimensions, and "writhing" points are also points, it is natural 
to ask whether there is any supporting evidence that the eigen- 
mode density is especially concentrated in point-like regions, 
rather than along lines, or surfaces, or 3-manifolds. 

A useful measure to quantify the localization of eigen- 
modes is the inverse participation ratio (IPR). The IPR of a 
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FIG. 7: Vortex correlation C^(N V ) for asqtad staggered eigenmodes on a 20 4 lattice at p/,iy = 3.3, a) full and b) center-projected configurations. 
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FIG. 8: Vortex correlation Cx(N v ) for overlap eigenmodes on a 16 4 lattice at fiiw = 3.3, a) full and b) center-projected configurations. 



normalized field pi(x) is defined as 

7 = iV£p?(*) 



(4.1) 



where N is the number of lattice sites x. Here, Pi(x) = (jc) 
and is the /-th, normalized (Y,xPi( x ) = 1)' lowest eigen- 
vector of the Dirac operator. The scaling of the IPR with lat- 
tice spacing is sometimes used to determine the dimension- 
ality of eigenmode concentration. Dimensionality is deduced 
from the IPR by reasoning that if the eigenmode has support 
mainly on a submanifold of dimension d, with a thickness in 
the 4 — d orthogonal directions which is a fixed number of 
lattice units, then the IPR should scale with lattice spacing 
as \/a A - d . This reasoning can lead to incorrect conclusions, 
because it is not necessarily true that the thickness of the lo- 
calization region is a constant number of lattice spacings, re- 
gardless of coupling. 

An instructive example is provided by the lowest eigen- 
mode of the covariant Laplacian in the adjoint representation, 



which was studied in Ref. [18]. In that case it was found that 
the IPR scaled like 1 /a 2 , suggesting an eigenmode concentra- 
tion on surfaces. Instead, it turns out that the lowest eigen- 
mode is sharply concentrated in a pointlike region. The pe- 
culiar scaling of the IPR arises because the volume b of the 
region of concentration, in lattice units, scales in a peculiar 
way. If this volume were a constant in physical units, then 
ba 4 would be constant. If instead the volume were constant 
in lattice units, then b itself would be constant. In fact, it is 
ba 2 which is constant; the volume of the eigenmode concen- 
tration region goes to zero in physical units, but infinity in lat- 
tice units, in the continuum limit. The naive deduction of the 
dimensionality of the concentration region, purely from the 
scaling of the IPR, leads in this case to an incorrect conclu- 
sion. For Dirac eigenmodes, conclusions based on the scaling 
of the IPR have not been entirely consistent with one another 
(cf. the overview in ref. [19]). Results of the MILC collabo- 
ration, with asqtad fermions, indicate a dimensionality d = 3 
[20], while the ITEP group has reported results, for overlap 
fermions and the Wilson action, consistent with d = [11]. 
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FIG. 9: Maximum density peak (center) of the first asqtad eigenmode on a 20 4 -lattice at = 3.3 with upper (above) and lower (below) 
z-slices of the same t-slice. Eigenmodes are computed on (a) full lattices, and (b) center-projected lattices. 



A third study, using overlap fermions and the Liischer-Weisz 
action, again indicates d = 3 [21] while the latest study of this 
group, using generalized IPR's defined using higher powers 
of pi(x), suggest eigenmode concentrations on manifolds of 
dimension between d = and d = 1 [22]. 

A possibly more reliable (if less quantitative) approach is 
to simply look at sample plots of p\(x) throughout the lat- 
tice volume. What we find, for the eigenmode density of the 
lowest eigenmodes of the asqtad and overlap operators, is that 
the eigenmode density is concentrated in very sharp peaks, in 
point-like regions of the lattice volume. 5 In Fig. 9 we display 
our data for Px(x), for the lowest eigenmode of the asqtad 
Dirac operator, in some two-dimensional slices of the four- 
dimensional lattice volume taken in the neighborhood of the 
point where px(x) is largest. In Fig. 9(a) we show the density 
of the lowest eigenmode computed from a typical configura- 
tion on the unprojected lattice, and in Fig. 9(b) we show the 
corresponding data for an eigenmode computed from a typi- 



cal center-projected lattice. Each lattice contains several sharp 
peaks of this kind; it is obvious that the concentration of eigen- 
mode density is in a point-like region, rather than being spread 
over a submanifold of higher dimensionality. In the figures we 
display a set of xy-plots of p(x,y,z,t) at various values of z and 
fixed ?, but there is an equally strong falloff of the peak as we 
move away from the maximum in the time direction. 

Figure 10(a) shows the same type of data for a zero mode of 
the overlap Dirac operator on 16 4 lattices, again evaluated on 
unprojected lattices generated from the Liischer-Weisz action 
at Plw = 3.3. Here again we find a handful of sharp peaks 
in the eigenmode density for any thermalized lattice configu- 
ration; one such peak is displayed in the plot, and it is con- 
centrated in a point-like region. The situation is very different 
for eigenmodes of the overlap operator evaluated in center- 
projected configurations. Instead of having a sharp peak, the 
eigenmode concentration in this case is very broad, extend- 
ing over most of the lattice volume, as seen in Fig. 10(b). 6 



5 The IPR itself is volume independent at fixed lattice spacing, indicating that 
the number of such peaks is proportional to the lattice volume, as expected. 



6 Note also the considerable difference in vertical scales, between Figs. 10(a) 
and 10(b). 
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FIG. 10: Maximum density peak (center) of the first overlap eigenmode on a 16 4 -lattice at p^iy = 3.3 with upper (above) and lower (below) 
z-slices of the same t-slice. Eigenmodes are computed on (a) full lattices, and (b) center-projected lattices (notice different scales!). 



However, we have already seen that the overlap operator, eval- 
uated on center-projected configurations, does not have any 
low-lying eigenmodes, let alone a zero mode. It is therefore 
not surprising that the eigenmode density is qualitatively dif- 
ferent from what is found in both the true zero modes of the 
overlap operator, and the "would-be" zero modes of the asqtad 
operator. 

Inspection of Figs. 9 and 10 indicate that eigenmode peaks 
are by far the sharpest for eigenmodes of the asqtad operator 
on the center-projected lattice. Taking account of the vertical 
scales in these figures, the peak in the asqtad-center-projected 
case (Fig. 9(b)), is about an order of magnitude higher than 
the peak in the asqtad-unprojected case (Fig. 9(a)). This dif- 
ference is of course reflected in a comparison of the IPRs of 
asqtad eigenmodes on the full and center-projected lattices, 
shown in Fig. 11, which indicate a far higher degree of eigen- 
mode concentration in the center-projected case. 

It is worth noting that on both the unprojected and center- 
projected lattices, the IPR rises as the lattice spacing becomes 
smaller, although the rate of increase is seen to be quite dif- 
ferent. According to Ref. [20] the IPR should go like 1 /a, in 
which case it should roughly double in going from $lw = 2.9 
to 3.3, Pliv = 3.1 to 3.5 or p^vv = 3.3 to 3.7. That behav- 



ior is at least crudely compatible with our data (Fig. 11(a)) 
for the unprojected lattices. As we have already stressed, this 
scaling of the IPR does not necessarily mean that eigenmode 
densities concentrate on three-volumes, and inspection of the 
density has revealed peaks in the eigenmode density in point- 
like regions. For center-projected configurations (Fig. 11(b)) 
the IPR of the lowest modes is roughly 1 1 at $lw = 2.9, 33 at 
Plw = 3.1, 140 at p LW = 3.3, and 400 at p LW = 3.5. Now, if 
the eigenmode density has support on pointlike regions whose 
volume is a fixed number of lattice spacings, regardless of 
Pliv, then the IPR should go like 1 /a 4 . That means the IPR 
should increase by about a factor of 4 from fiiw = 2.9 to 3.1, 
from p/,iv = 3.1 to 3.3 and from pny = 3.3 to 3.5. This is not 
so far off the actual results. 

It is not too surprising that would-be zero modes of the asq- 
tad operator would be very highly concentrated when evalu- 
ated on center-projected lattices. On thin vortices, topological 
charge is concentrated not just in point-like regions, but in fact 
at individual lattice sites on the dual lattice, where thin vortex 
sheets writhe and/or intersect. Since zero modes concentrate 
on regions of non-zero topological charge density, and topo- 
logical charge is concentrated at individual sites on the center- 
projected lattice, the high degree of localization of the lowest- 
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(a) full (unprojected) lattices (b) center-projected lattices 

FIG. 11: Inverse Participation Ratio vs. eigenmode number for asqtad eigenmodes, computed on a 20 4 -lattice at a variety of $lw values for: 
a) unprojected lattices, and b) center-projected lattices. 



lying modes, in volumes of lattice-scale extension, is to be 
expected. On unprojected lattices the sources of topological 
charge, whether vortices, instantons, calorons, or something 
else, are more spread out, and there is no particular reason 
to expect that the eigenmode density would concentrate in 
tiny regions of lattice scale extension. What we would like to 
know, of course, is which of the candidate sources of topolog- 
ical charge density is giving the main effect. Our density plots 
give a strong indication that the charge density concentrates in 
point-like regions, rather than surfaces or three-volumes, but 
this fact would be compatible with instanton, caloron, and vor- 
tex (intersection/writhing) sources. The eigenmode density- 
vortex surface correlations, seen in Figs. 7 and 8 provide a 
modest degree of evidence in support of a vortex origin of 
topological charge density. There is also the possibility, of 
course, that topological charge density may come from more 
than one type of source. 



V. CONCLUSIONS 

We have shown that center-projected SU (2) lattice config- 
urations give rise to a dense low-lying Dirac eigenvalue spec- 
trum, as required for chiral symmetry breaking, for a massless 
lattice Dirac operator in the asqtad formulation. In contrast, 
this low-lying spectrum is not found for chirally improved 
and overlap Dirac operators on center-projected lattices, for 
reasons which are almost certainly connected to the lack of 
smoothness of center-projected configurations. Chiral sym- 
metry is absent in the chirally-improved Dirac operator on 
such configurations, while for the overlap operator an exact 
symmetry is present, but is strongly field-dependent for rough 
configurations (and thus quite different from the continuum 
symmetry). In the case of the overlap operator we have found 
that a moderate degree of "smoothing" of the center-projected 
lattice brings back the low-lying spectrum. In a staggered for- 
mulation such as asqtad, the smoothness of the lattice config- 



uration has nothing to do with the exact chiral symmetry, and 
the low-lying modes are present for thin vortex configurations. 

There is a general expectation, based on the old Casher ar- 
gument [2], that gauge field configurations with the confine- 
ment property ought to also break chiral symmetry. Our re- 
sults indicate that this expectation holds for confining ensem- 
bles of thin center vortices, at least when the relevant lattice 
Dirac operator (asqtad) has some subset of the exact chiral 
symmetry required by the Casher argument. Although non- 
confining lattice configurations may also break chiral sym- 
metry, we find that vortex-removed configurations (which are 
non-confining) have a large gap in the asqtad Dirac spectrum 
around A, = 0, indicating an unbroken chiral symmetry. 

We have also looked at the correlation of vortex location 
with the densities of low-lying Dirac eigenmodes, following 
the earlier work of Kovalenko et al. [13]. We find that for both 
asqtad and overlap eigenmodes computed on an unmodified 
lattice, there is a significant positive correlation between low- 
lying eigenmode densities and the location of thin vortices on 
the corresponding center-projected lattice. It is found that this 
correlation is greatest in the neighborhood of points where a 
large number of vortex plaquettes meet, such as would be the 
case for vortex intersections, "writhings", or any combination 
of these effects. Since thin vortices in D = 4 dimensions in- 
tersect and writhe only at points, we would therefore expect 
that low-lying Dirac eigenmodes are especially concentrated 
in point-like regions. A simple inspection of eigenmode den- 
sity distributions reveals that these densities do indeed possess 
sharp peaks in point-like regions. The vortex-eigenmode cor- 
relations, and the peaks in the eigenmode densities at point- 
like regions, together provide a degree of support for the pic- 
ture advanced by Engelhardt and Reinhardt [14], in which 
topological charge would tend to concentrate in the neighbor- 
hood of center vortex intersection and writhing points. 

It is of interest to compare our results to some related work 
by Ilgenfritz et al. [22]. These authors report that the high 
density regions of overlap zero modes are concentrated in re- 
gions of dimensionality in the range of zero to one; the lower 
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bound is consistent with our finding (for both overlap and asq- 
tad modes) of sharp peaks in the eigenmode density, located 
in point-like regions. We also take note of the earlier work 
of Gubarev et al. in ref. [11], which found IPR's for over- 
lap eigenmodes consistent with concentration in point-like re- 
gions (for an overview of this and other IPR results, cf. [19]). 
Ilgenfritz et al. [22] also find a significant positive correlation 
between the topological charge density at a lattice site, and 
the probability that a given site is adjacent to a P-vortex sheet 
or monopole line on the dual lattice. That result is consistent, 
assuming a correlation of zero-mode density and topological 
charge density, with the correlations that we (and Kovalenko 
et al. [13]) have found between vortices and the densities of 
low-lying Dirac eigenmodes. 

To summarize: there are significant correlations between 
center vortices and the low-lying modes of both the asqtad 
and overlap Dirac operators, and this correlation steadily in- 
creases with vortex connectivity (N v ). We also find that the 
thin vortices found in center projection give rise to a low-lying 
spectrum of Dirac eigenmodes, as required for %SB by the 
Banks-Casher formula, providing that the chiral symmetry of 
the Dirac operator does not depend on the smoothness of the 
lattice configuration. Vortex-removed configurations do not 
have these low-lying eigenmodes, and therefore do not break 
chiral symmetry. Taken together, these results indicate that 
center vortices have a strong effect on the existence and prop- 
erties of low-lying eigenmodes of the Dirac operator. 



APPENDIX 



coefficients c, = + 47COCoA; for one-loop corrections A, have 
been computed by Liischer and Weisz for both St/ (2) and 
SU(3) (Table 1 in Ref. [23]). The continuum limit behav- 
ior of the Liischer- Weisz action can be further improved by 
making the lattice links more "continuum like". At the mean 
field level this entails setting E/„ — > Un, where one possible 
choice for the mean field (or "tadpole") factor uq is using the 
expectation value of the average plaquette 



--< -ReTrf/„, >^ 4 
N 1 



(A.2) 



The Liischer- Weisz action can now be tadpole improved by 
explicitly pulling a Uq 1 factor out of each link and replacing 
oco in the one-loop perturbatively renormalized coefficients 
c, with a nonperturbatively renormalized coupling oc 4 defined 
through [24] 



-4^ with 



: 0.36626271 



N 2 -l 
N 



1.72597, for Af = 2 
3.06839, forAf = 3 



(A.3) 



Defining $lw = u Q 4 $ (since U p i involves 4 links) the im- 
proved action reads for SU (2) [24] 



S = frwVSpi - ^[1 + 0.2227%] £S,t 



pi 



0.02224 ^a*£S P g 



(A.4) 



The gauge action used in this work is a tadpole improved 
version of the one-loop continuum limit improved SU (2) ac- 
tion of Liischer and Weisz [23, 24]. 



a) 



b) 



O 



-A 



FIG. 12: Liischer- Weisz action Wilson loops: a) standard plaquette, 
b) 2 x 1 rectangle and c) 1 x 1 x 1 parallelogram 

The standard Liischer- Weisz action removes leading ultra- 
violet cutoff effects by adding a few next-to-nearest neighbour 
terms to Wilson's action. In addition to the standard plaquette 
(labeled "pi") term, it includes a sum over all 2 x 1 (planar) 
rectangle (labeled "rt") and over all 1 x 1 x 1 parallelogram 
(labeled "pg") Wilson loop terms (see Fig. 12). For SU (N) 
lattice gauge fields U^(x) living on a four-dimensional (p = 
0, 1,2,3) hypercubic lattice with sites x and lattice spacing a, 
the improved action reads 



X y pi rt pg 



(A.l) 



where p denotes the (inverse) coupling constant and Si = 
jjRe. Tr(l — U) with [/, the corresponding Wilson loops. The 



The tadpole factor mo is determined during thermalization and 
then kept fixed. In view of extracting physical quantities, 
we fit the time-dependent potential V(R,T) = log(W[R,T — 
\}/W[RJ]) with W[R,T] the Wilson loop of size R x T 
in space-/time-direction respectively at some fixed T, to an 
ansatz V(R) = OjR — c/R + vq (linear-plus-Coulomb fit). In 
order to obtain an asymptotic lattice string tension a/ ur we fit 
the extracted string tensions Ot for several T values to some 
stabilizing function f(T) = e\p(—kT + d) + a with a giving 
the asymptotic (T — > °°) value. All fits were done by least- 
square routines. To set the scale we use the physical string 
t/a = ^/Gp~g~ s w 0.44GeV [25], to determine the 



tension, ^/Cu„/a = ^/a^J 
lattice spacing a. Table I lists the data for runs of string ten- 
sion determination on 20 4 -lattices with a 1000 thermalization 
steps, 1000 measurements separated by 200 iterations each. 



pLW 


a lot 


a[fm] 


2.9 


0.3756 ±0.0053 


0.2749 ±0.0019 


3.1 


0.2254 ±0.0033 


0.2129 ±0.0016 


3.3 


0.1112 ±0.0017 


0.1495 ±0.0012 


3.5 


0.0635 ±0.0007 


0.1138 ±0.0006 


3.7 


0.0401 ±0.0003 


0.0898 ±0.0003 


4.0 


0.0225 ±0.0002 


0.0673 ±0.0002 



TABLE I: Lattice string tension ai at and lattice spacing a, as ex- 
tracted from linear-plus-Coulomb fits to V(R, T = fixed) 
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(a) P- vortex surface density vs. coupling constant puy- "Two loop" line is 
the scaling prediction with y/p/ 6A 2 = 50. 
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(b) Center-projected Creutz-ratios at poy =3.1 — 3.5. Horizontal bands 
indicate asymptotic string tensions on the unprojected lattice. 



FIG. 13: Evidence for: a) asymptotic scaling of vortex densities; and b) center dominance for the string tension. 
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(b) Vortex removal. 



FIG. 14: Effects of vortex restriction, and vortex removal, on Wilson loops, (a) Vortex limited Wilson loops W n on 20 4 lattice at fiiw = 3.1. 
Note the approach, as loop area increases, to the limit (— l)"W„. (b) The effect of vortex removal. 20 4 lattice at $lw = 3.3 for full, center- 
projected and vortex-removed data. Horizontal bands indicate asymptotic string tension on the unprojected lattice with errorbars. 



All previous numerical checks of the vortex mechanism, 
via maximal center gauge fixing and center projection, have 
been carried out on thermalized lattices generated from Monte 
Carlo simulations of the Wilson action (a discussion of these 
tests and their significance can be found in Ref . [ 1 ] .) It is im- 
portant to repeat these checks for the new action. It turns out 
that there are no surprises, and vortex results derived using the 
Luscher-Weisz action are consistent with the previous work 
in the Wilson action. Here the gauge-fixing procedure is the 
over-relaxation method applied to one configuration, resulting 
in a single (random) gauge copy. No attempt is made to find 
the true global maximum of squared link variables, or the best 
out of a set, so technically the procedure corresponds to av- 
eraging over all gauge copies of direct maximal center gauge 
in the Gribov region. We have compared this procedure to the 
"best copy of five" prescription, and have found no significant 
difference in the results. 



Figure 13(a) shows the P-vortex density (in lattice units) vs. 
Pz,w- The solid line is the two-loop asymptotic freedom be- 
havior for this quantity, for a choice of vortex density p (area 
per unit volume) in physical units satisfying yp/ 6 A 2 = 50. 
Figure 13(b) is a test of vortex dominance; the figure shows 
the values for Creutz ratios %(R,R) on the center-projected lat- 
tice at various R and various couplings. The horizontal bands 
indicate the asymptotic string tension on unprojected lattices 
(middle lines), together with their corresponding error bars 
(upper and lower lines). The center-projected Creutz ratios 
approach the asymptotic string tension, but at largest available 
R values are still lower by about 10%. 

In Fig. 14(a) we show our data for the ratios of vortex- 
limited Wilson loops Wi/Wo and W2/W0. A vortex-limited 
Wilson loop is a loop built of unprojected links, which is 
pierced by n P-vortices on the corresponding center-projected 
lattices. It is expected that for sufficiently large loops in 
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SU(2) gauge theory, Wi/Wo ->• -1 and W 2 /W -»■ 1, at least 
if the vortex piercings are near the middle of the loop, rather 
than lying near the perimeter. The data shown in Fig. 14(a) 
was collected with this restriction on P-vortex piercings. 
Finally, the effect of vortex removal on Creutz ratios is shown 
in Fig. 14(b), together with full and center-projected results, 
at Plw = 3.3. All of this data is similar to corresponding 
results previously obtained from simulations of the SU (2) 
Wilson action. 
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